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ABSTRACT 

The source responsible for reionizing the universe at z > 6 remains uncertain. While an energetically 
adequate population of star-forming galaxies may be in place, it is unknown whether a large enough 
fraction of their ionizing radiation can escape into the intergalactic medium. Attempts to measure this 
escape-fraction in intensely star-forming galaxies at lower redshifts have largely yielded upper limits. 
In this paper we present new HST COS and archival FUSE far-UV spectroscopy of a sample of eleven 
Lyman Break Analogs (LB As), a rare population of local galaxies that strongly resemble the high-z 
Lyman Break galaxies. We combine these data with SDSS optical spectra and Spitzer photometry. 
We also analyse archival FUSE observations of fifteen typical UV-bright local starbursts. We find 
evidence of small covering factors for optically-thick neutral gas in three cases. This is based on two 
independent pieces of evidence: a significant residual intensity in the cores of the strongest interstellar 
absorption-lines tracing neutral gas and a small ratio of extinction-corrected Ha to UV plus far-IR 
luminosities. These objects represent three of the four LBAs that contain a young, very compact 
(~ 10^ pc), and highly massive (~ 10^ M©) dominant central object (DCO). These three objects also 
differ from the other galaxies in showing a significant amount of blueshifted Lya emission, which may 
be related to the low covering factor of neutral gas. All four LBAs with DCOs in our sample show 
extremely high velocity outfiows of interstellar gas, with line centroids blueshifted by about 700 km 
s~^ and maximum outflow velocities reaching at least 1500 km s^^. We show that these properties 
are consistent with an outflow driven by a powerful starburst that is exceptionally compact. We 
speculate that such extreme feedback may be required to enable the escape of ionizing radiation from 
star forming galaxies. 

Subject headings: galaxies: evolution — galaxies: intergalactic medium — galaxies: high redshift — 
galaxies: ISM — galaxies: kinematics and dynamics — galaxies: starburst 



1. INTRODUCTION 

Following the epoch of recombination at z ^1000, the 
next major event in the history of the universe was 
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its reionization by an early generation of massive stars 
and/or black holes. One of the major goals of modern 
cosmology is to understand this process: what sources 
were responsible and when did it occur? Observations 
of the cosmic microwave background with WMAP show 
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that reionization occurred at a mean redshift of 11. Oil. 4 
(Dunkley et al. 2009), while observations of absorption 
due to the Lyman a forest in quasars show that reion- 
ization was not quite complete until a redshift of about 
six (e.g. Fan et al. 2006). 

Based on the observed luminosity functions of active 
galactic nuclei and star-forming galaxies at the relevant 
redshifts, the latter appear to be the more energetically 
plausible source for reionization (e.g. Bouwens et al. 
2010). However, a major uncertainty in the contribu- 
tion of star forming galaxies is the poorly constrained 
fraction of ionizing photons that arc able to leak out of 
such galaxies into the intcrgalactic medium. Boiiwcns et 
al. (2010) estimate that an escape fraction of 20% (60%) 
is required for consistency with the WMAP results at the 
2(7 (1 ct) level. 

Many attempts have been made to mcasiirc this 
escape-fraction in star-forming galaxies at low-redshift 
(Lcithcrcr et al. 1995; Hurwitz et al. 1997; Dcharveng et 
al. 2001; Heckman et al. 2001 - hereafter HOI; Bergvall 
et al. 2006; Grimes et al. 2009 - hereafter G09), interme- 
diate redshift (Malkan et al. 2003, Siana et al. 2007,2010; 
Cowie et al. 2009) and high- redshift (Steidel et al. 2001; 
Giallongo et al. 2002; Fernandez- Soto et al. 2003; Inoue 
et al. 2005; Vanzella et al. 2010). For the most part, 
these measurements have only yielded upper limits, al- 
though significant escaping Lyman continuum emission 
has been detected in 17 out of a sample of 198 Lyman 
Break and Lyman a galaxies at z > 3 (Iwata et al. 2009), 
consistent with earlier results for a smaller sample of Ly- 
man Break galaxies (Shapley et al. 2006). In contrast, 
Vanzella et al. (2010) found evidence for escaping Ly- 
man continuum in only 1 of 100 Lyman Break galaxies 
at z ~ 4. 

Why do only a minority of the high-z galaxies and 
none of the low-z galaxies show a significant leakage of 
ionizing radiation? In this regard, it is instructive to 
note that an HI column density of only 1.6 xlO^^ cm~^ 
is required to produce unit optical dcipth at the Lyman 
edge. The area-averaged column densities of cold gas 
in star-forming galaxies are much larger, ranging from 
10^^ cm^^ in ordinary low-rcdshift galactic disks to as 
high as 10^^ cm~^ in local starbursts and high-z star- 
forming galaxies (Kennicutt 1998; Gcnzcl ct al. 2010). 
The escape of ionizing radiation must then be determined 
by the topology of the interstellar medium, which will 
in turn be strongly influenced by feedback from massive 
stars. It may be that rather exceptional circumstances 
are required to produce the required channels of very low 
column-density neutral gas. Such processes have been 
extensively studied from a theoretical perspective ovre 
the past decade (e.g. Dove et al. 2000; Clarke & Oey 
2002; Fujita et al. 2003; Gnedin & Kravtsov 2008; Wise 
& Cen 2009; Razoumov & Somer-Larson 2010; Yajima 
et al. 2010). 

These considerations provide the motivation of the 
present paper: to find local examples of galaxies in which 
a potentially significant fraction of the ionizing radiation 
may be escaping, and to then exploit the relative ease of 
observing such bright local objects to gain insight into 
the physical processes that facilitate this escape. To be 
able to assess the cosmological implications of the re- 
sults of such an approach, it is important to observe local 
galaxies that most closely resemble the type of UV-bright 



high-redshift galaxies that are candidates for the reion- 
ization of the universe. 

We began such an investigation in Overzier et al. (2009 
- hereafter O09), where we reported on Spitzer photome- 
try, HST imaging, and Sloan Digital Sky Survey (SDSS) 
optical spectroscopy of a sample of 31 Lyman Break 
Analogs (LB As). These LB As are a very rare popula- 
tion of local galaxies that strongly resemble the high-z 
Lyman Break galaxies (Heckman et al. 2005; Hoopes et 
al. 2007; Overzier at al. 2008;2010; Basu- Zych et al. 
2007;2009). They are defined to have a range in far-UV 
luminosity (> 2 x IQ^^ Lq) and far-UV effective surface 
brightness (> 10^ Lq kpc~^) that are similar to typical 
Lyman Break galaxies. 

We showed in O09 that the morphologies of LBAs as 
revealed in multi-band HST images fall into two distinct 
categories. The first class (25 objects) consists of galaxies 
in which multiple bright complexes of clumps and knots 
reside in a galaxy whose disturbed morphology is sugges- 
tive of a major merger or strong tidal interaction. The 
second class (six objects) consists of galaxies whose UV 
image is dominated by a single extremely massive (Af* ~ 
one-to-several billion Mq) and compact (radius of-order 
100 pc) central object (a dominant central object, or 
DCO). We pointed out that the star- formation rates de- 
rived for these six galaxies based on the extinction- cor- 
rected Ha luminosities were systematically smaller than 
those derived from the Spitzer mid-IR luminosity (by 
factors of ~ 2 to 5), and speculated that a possible ex- 
planation was that these galaxies were leaking ionizing 
radiation into the intergalactic medium. 

In this paper we test this idea using new Hubble Space 
Telescope (HST) far-UV observations with the Cosmic 
Origins Spectrograph (COS) of a sample of eight LBAs 
from O09. We will combine these new data with archival 
Far Ultraviolet Spectroscopic Explorer (FUSE) data for 
three other LBAs and for fifteen more typical UV-bright 
local starbursts. We will describe the sample selection 
in section 2 and the observations and data analysis in 
section 3. We will present our results and consider their 
interpretation in section 4, and summarize these results 
and their implications in section 5. 

2. SAMPLE SELECTION 

We have previously used the Hopkins Ultraviolet Tele- 
scope and FUSE to search for the direct evidence of the 
escape of ionizing radiation from local starbursts (Lei- 
therer et al. 1995; HOI; G09), but the modest sensitiv- 
ity of these facilities was sufficicmt only to observe three 
LBAs (none of which contained a DCO). With the recent 
installation of the Cosmic Origins Spectrograph (COS) 
on the Hubble Space Telescope (HST), we are now in a 
position to obtain high-quality far-UV spectra of a sig- 
nificantly larger sample of these galaxies, including some 
with DCOs. 

For the redshifts of the LBAs in our HST program 

[z ^ 0.1 to 0.25), the sensitivity of COS is poor in the 
rest-frame Lyman continuum. The sensitivity is much 
better at longer wavelengths, so in this paper we follow 
HOI and G09 and use the residual intensity in the core 
of the strongest far-UV interstellar absorption-lines from 
the neutral phase of the interstc^Uar medium as a probe 
of the intensity in the region below the Lyman edge. We 
will describe this in more detail in section 3 below. 
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Fig. 1. — HST rest-frame optical images of the eight Lyman Break Anologs observed with COS. The circle indicates the 2.5 arcsec 
aperture of COS. For details on the images see Overzier et al. (2009). Note the diffraction spikes visible in the images of 0213+12, 
0808+39, 0921+45, and 2103-07 produced by the barely-resolved Dominant Central Object. 



This approach requires the observations to be obtained 
through an aperture that encloses the bulk of the star- 
burst's emission (so that a global constraint is obtained). 
It also requires that the width of instrumental spectral 
line-spread function is significantly smaller than the char- 
acteristic velocity dispersion in the starburst (so that the 
interstellar absorption lines are well resolved). Finally, 
the data need to have a high signal-to-noise ratio in the 
far- UV continuum so that useful constraints are derived. 

Based on these considerations, we will consider two 
data sets in this paper. The first consists of the sample 
of eighteen galaxies studied by HOI and G09. These are 
local starburst or star-forming galaxies observed with the 
Far-Ultraviolet Spectroscopic Explorer (FUSE - Moos et 
al. 2000), each having signal-to-noise better than 4.6 per 
0.078 A (~ 25 km s'^) spectral element in the FUSE 
LiFlA channel (see G09). The observations and prop- 
erties of this sample are described in Table 1. Three of 
these galaxies are in fact LBAs without a DCO: Mrk 54 
(Deharveng et al. 2001), Haro 11 (Grimes et al. 2007), 
and VV 114 (Grimes et al. 2006). However, most of 
these galaxies have considerably lower UV luminosities 
and star-formation rates than the LBAs. They also span 
a much broader range in galaxy mass and metallicity (see 
G09 and O09, and compare Table 1 and 2). 

The second is a sample of eight LBAs observed with the 
Cosmic Origins Spectrograph (Froning & Green 2009) 
on the Hubble Space Telescope (Program 11727: PI T. 
Heckman). These are members of a sample of 31 LBAs 
with HST UV images discussed by O09 and were selected 
for spectroscopy based on a high UV flux through the 
COS aperture and a compact UV size (so that the COS 
line-spread function is not significantly degraded). These 
observations and the properties of this sample are listed 
in Table 2 and HST images are shown in Figure 1. In the 
discussion to follow we will refer to the eight LBAs with 



HST COS data and the three LBAs with FUSE data 
as the LBA sample. We will refer to the other fifteen 
galaxies with FUSE data as the local starburst sample. 

3. OBSERVATIONS AND DATA ANALYSIS 

3.L The Data 

The observations and data reduction of the FUSE sam- 
ple have been described in detail in HOI and G09. We 
refer the reader to these papers. All the data were ob- 
tained through the 30 x 30 arcsec LWRS aperture, ex- 
cept for the cases of NGC 5253 and NGC 7714 (which 
used the 4 x 20 arcsec MDRS aperture). As shown by 
G09 these apertures encompass most or all of the star- 
burst in the far-UV. These spectra cover the observed 
wavelength range from 905 to 1187 A ( see Table 1 for 
the corresponding range in the rest-frame). Depending 
on the angular size of the starburst in the far-UV, the 
instrumental spectral resolution is i? 5000 to 14,000, 
corresponding to a velocity dispersion of cr ~ 9 to 25 km 
s~^ (G09). In all cases the interstellar absorption lines 
in the starbursts are well resolved. 

For the HST-COS sample we have used the COS 
G130M and G160M gratings to obtain spectra of our 
eight targets. As can be seen in Figure 1, the COS aper- 
ture encompasses most or all of the galaxy. We have 
retrieved these data from the HST MAST archive af- 
ter they have been processed through the standard COS 
pipeline. The merged spectra cover a range from about 
1160 to 1780 A with the corresponding range in the rest- 
frame wavelength given for each sample member in Ta- 
ble 2. As with the FUSE data, the instrumental spectral 
resolution depends on the angular size of the far-UV con- 
tinuum within the COS aperture. Our HST UV images 
(O09) yield sizes of 0.1 to 0.5 arcsec (FWHM), and we es- 
timate that the resulting spectral resolution is i? ^ 3000 
to 13,000, corresponding to a velocity dispersion of u ~ 
10 to 43 km s~i (similar to the FUSE data). 




Fig. 2.— The full COS G130M and G160M spectrum of LBA0921+45 is shown in black, while an illustrative Starburst99 model spectrum 
is shown in red (constant star formation for ton Myr with 0.4 solar metallicity and a Chabrier IMF). Both spectra have been normalized 
to unit flux (see text). Prominent spectral features in the starburst arc labeled at the top, while Milky Way and geocoronal features are 
labeled at the bottom. 



For the new COS spectra we have fit each merged 
G130M and G160M spectrum with a Starburst99 (Lei- 
therer et al. 1999) model spectrum for a young stel- 
lar population with a metallcity similar to that of the 
corresponding LB A (Table 2). We did not attempt to 
obtain a rigorous best-fit, but simply adopted a stan- 
dard Chabrier (2003) IMF and assumed a constant rate 
of star-formation. We varied the duration of the star- 
formation until we obtained (by eye) a satisfactory fit to 
the strong stellar wind features. We then took the ratio 
of the model and the data and fit the result with a low- 
order polynomial (avoiding spectral regions with strong 
emission or absorption features). We then divided the 
data by the polynomial to obtain a normalized spectrum. 
An example of this is shown in Figure 2. 

3.2. The Analysis Approach 

The opacity of the interstellar medium to Lyman con- 
tinuum photons is due to two primary sources: the pho- 
toelectric opacity of the neutral hydrogen and the opac- 
ity of dust. The latter significantly affects the entire UV 
continuum in typical local starbursts and high-z galax- 
ies. We will focus first on measuring constraints on the 
photoelectric opacity. Later we will comment on the im- 
portance of dust. 

The only direct probe of the photoelectric opacity is 
the direct measurement of the relative intensity in the 
continuum shortward of the Lyman edge. The FUSE 
data provide clean access to this spectral region and ad- 
equate signal-to-noise in only 5 targets. G09 have used 
these data to derive upper limits to the escape of ion- 
izing radiation for these galaxies (see also Deharveng et 
al. 2001; Bergvah et al. 2006). The COS data do not 
sample this spectral region in any galaxy. Thus, to gain 
insight into the escape of ionizing radiation for our two 
samples we will follow HOI and G09 and use as probes 
the residual intensity in the core of the strongest inter- 
stellar absorption lines that primarily arise in the neutral 
interstellar medium. 

For the FUSE data the strongest such line is 



CIIA1036.3 and for the HST-COS data the strongest such 
line is CIIA1334.5. Thus in both samples we will use the 
same ionic species as a tracer of the neutral ISM. ^ We 
have measured the relative residual intensity in the core 
of the CII absorption line, following the same methodol- 
ogy as in HOI and G09. The results are listed in Tables 
1 (FUSE) and 2 (COS). 

4. RESULTS & INTERPRETATION 

4.1. The Ionization Source 

Before considering constraints on the escape of ionizing 
radiation, it is important to establish that an ionizing 
population of massive stars (O stars) is in fact present 
in our targets. The strongest spectroscopic signature of 
these stars is provided by the broad P-Cygni profiles in 
the resonance transitions of highly ionized species that 
arise in the stellar winds (e.g. Robert et al. 1993; 2003). 

The strongest such feature that is present in all our 
COS spectra is NVA1240. In Figure 3 we show the NV 
profile for our eight targets and overplot the generic ex- 
amples of Starburst99 model spectra described above. 
We defer to a future paper a detailed discussion of the 
properties of the stellar populations in the LBAs. Here 
we merely point out that the far-UV spectra of the LBAs 
are all consistent with a normal population of ionizing 
stars expected for a strong starburst. 

4.2. Constraints on the Escape of Ionizing Radiation 

In Figure 4 we show the CIIA1334.5 absorption-line 
profiles for our COS data on eight LBAs. The corre- 
sponding profiles for the CIIA1036.3 line from the FUSE 

^ In principle, a significant fraction of the gas- phase carbon 
in the neutral medium could be in the form of CI (which has an 
ionization potential of 11.52 eV). The CIA1277.2 line should be 
comparable in strength to the two CII lines in this case, however it 
is extremely weak in our COS spectra and in the spectra of other 
local starburst galaxies (e.g. Heckman et al. 1998; Leitherer et 
al. 2002; Schwartz et al. 2006; Leitherer et al. 2010) and high-z 
Lyman Break galaxies (Shapley et al. 2003; Pettini et al. 2002; 
Steidel et al. 2010; Quider et al. 2009;2010; Dessauges-Zavadsky 
et al. 2010). 
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TABLE 1 
FUSE Sample 



Name 


LBA 


^sys 


Range 


^ esc,rel 


^ esc,abs 


SFR 


log M« 








Y 


6175 


887 


to 


1163 


< 0.08 


< 0.006 


26 


10.2 


-180 


-450 


VV114 


Y 


6016 


887 


to 


1164 


< 0.02 


<0.001 


66 


10.8 


-150 


-650 


NOr 1 1 do 


N 


1501 


900 


to 


1181 


< 0.04 


< 0.014 


0.83 


9.4 


-30 


-400 


ooo uooo uo^ 


N 


4043 


893 


to 


1171 


< 0.18 


< 0.18 


0.32 


7.8 


30 


-130 


TOT, 0440-^81 


N 


12291 


869 


to 


1140 


< 0.04 


< 0.016 


5.0 


10.0 


-120 


-450 


NOP 1 7n^ 


N 


632 


903 


to 


1184 


< 0.01 


< 0.008 


0.16 


8.6 


-12 




NOC 1 74-1 


N 


4037 


893 


to 


1171 


< 0.01 


< 0.002 


6.0 


9.7 


-35 


-500 


IZW 18 


N 


751 


903 


to 


1184 


< 0.1 


< 0.1 


0.016 


7.1 


-3 


-150 


NGC 3310 


N 


983 


902 


to 


1183 


< 0.05 


< 0.005 


2.8 


9.8 


-170 


-600 


HARO 3 


N 


963 


902 


to 


1183 


< 0.005 


< 0.001 


0.41 


8.9 


10 


-150 


NGC 3690 


N 


3121 


896 


to 


1175 


< 0.01 


< 0.001 


82 


10.9 


-120 


-900 


NGC 4214 


N 


291 


904 


to 


1186 


< 0.01 


< 0.002 


0.13 


9.0 


-40 


-250 


MRK 54 


Y 


13450 


866 


to 


1136 


< 0.09 


< 0.05 


21 


10.4 


-90 


-550 


M 83 


N 


517 


903 


to 


1185 


< 0.06 


< 0.01 


3.5 


10.7 


9 




NGC 5253 


N 


405 


904 


to 


1185 


< 0.06 


< 0.02 


0.33 


9.3 


17 




IRAS 19245-4140 


N 


2832 


897 


to 


1176 


< 0.21 


<0.11 


2.1 


9.1 


-60 


-350 


NGC 7673 


N 


3392 


895 


to 


1174 


< 0.1 


< 0.017 


4.8 


10.0 


-80 


-350 


NGC 7714 


N 


2803 


897 


to 


1176 


< 0.03 


< 0.01 


6.9 


10.2 


-60 


-500 



Note. — Column 1: Galaxy name. Column 2: Is the galaxy a Lyman Break Analog? See O09 for a definition. 
Column 3: Galaxy systemic velocity (km s~^) from G09 or HOI. Column 4: Wavelength range covered by the 
FUSE data in the galaxy rcst-framc (A). Column 5: The relative escape fraction of ionizing radiation based on 
the residual relative intensity in the CIIA1036.3 line - sec text for details. Column 6: The absolute escape fraction 
including the effect of attenuation due to dust opacity - sec text for details. Column 7: The star formation rate 
(M.0 year""'"), based on the sum of the total infrared and far-ultraviolet luminosities and assuming a Chabricr 
2003 IMF - sec GOO for details. Column 8: The log of the galaxy stellar mass (Mq) based on the 2MASS K-band 
luminosity and assuming a K-band mass-to-light ratio of 0.4 (Bell & de Jong 2001). Column 9: The centroid 
velocity of the strong interstellar absorption- lines with respect to the galaxy systemic velocity (km sec~^). See 
G09 for details. Column 10: The maximum outflow velocity (in km sec~^), following Steidel et al. (2010). 



TABLE 2 
COS Sample 



Name 


DCO 


'Z'sys 


Range 




esc,abs 


SFR/ij+f/v 


SFRhc 


log M, 








/3 


[O/H] 


0055-00 


N 


0.16744 


970 to 1513 


<0.05 


<0.01 


29 


23 


9.7 


-230 


-640 


-0.21 


-1.5 


-0.41 


150+13 


N 


0.14668 


1022 to 1572 


<0.1 


<0.01 


56 


19 


10.3 


-150 


-480 


-0.78 


-1.6 


-0.30 


0213-1-12 


Y 


0.21902 


945 to 1440 


0.35 


0.05 


55 


6.7 


10.5 


-770 


-1500 


0.36 


-0.9 


0.05 


0808-1-39 


Y 


0.09123 


1038 to 1607 


0.3 


0.12 


11 


3.7 


9.8 


-680 


-1500 


0.27 


-1.0 


0.05 


0921+45 


Y 


0.23499 


948 to 1430 


0.3 


0.04 


84 


26 


10.8 


-650 


-1.500 


0.75 


-1.4 


-0.02 


0926+45 


N 


0.18072 


968 to 1515 


0.6: 


0.4; 


10 


13 


9.1 


-300 


-700 


0.14 


-2.3 


-0.60 


0938+54 


N 


0.10208 


1029 to 1614 


<0.1 


<0.05 


13 


13 


9.4 


-90 


-610 


-0.17 


-2.1 


-0.50 


2103-07 


Y 


0.13689 


1030 to 1585 


<0.05 


<0.01 


63 


41 


10.9 


-720 


-1500 


0.04 


-0.5 


0.01 



Note. — Column 1: Galaxy name; Column 2: Docs the galaxy contain a Dominant Central Object (DCO)? See O09 for details. Column 3: The 
galaxy redshift from Overzier et al. (2011). Column 4: The wavelength range covered by the COS data in the galaxy rest-frame (A). Column 5: The 
relative escape fraction of ionizing radiation based on the residual relative intensity in the CIIA1334.5 line. In the case of 0926+45 the CII line has a 
significant residual intensity, but there is no deficiency of Ha emission (unlike the cases of 0213+12, 0808+39, and 0921+45). Only these three latter 
objects are good candidates for the escape of ionizing radiation. See text for details. Column 6: The absolute escape fraction including the effect 
of attenuation due to dust opacity - see text for details. Column 7; The star formation rate (Mq year~^), based on the sum of the total infrared 
and far-ultraviolet luminosities and assuming a Chabrier 2003 IMF - see Overzier et al. 2011 for details. Column 8: The star formation rate {Mq 
year~^), based on the extinction-corrected Ha luminosity measured from the SDSS spectra - see Overzier et al. 2011 for details. Column 9; The log 
of the galaxy stellar mass (Mq). Sec O09 for details. Column 10; The centroid velocity of the SiIIIA1206.5 interstellar absorption-line with respect 
to the galaxy systemic velocity (km scc^^). Column 11: The maximum outflow velocity seen in the SiIIIA1206.5 line (in km sec^^), following Steidel 
et al. (2010). Column 12: The ratio of the equivalent widths of the blucshifted and rcdshil'tcd components of the Lya line. Values near one imply 
significant blucshifted Lya emission. See text for details. Column 13: The spectral slope of the UV continuum defined by Fx oc A^. An unrcddcned 
starburst will have /3 ~ -2.5 to -2. Six of the eight LBAs suffer significant dust reddening. Column 14: The log of gas-phase O/H ratio relative to 
solar (Asplund et al. 2009). See Overzier et al. (2009) for details. 

data can be found in HOI and G09. A striking re- 
sult is that there is a significant residual intensity in 
the absorption-line profiles in four LBAs (LBA0213-I-12, 
LBA0808-h39, LBA0921-h45, and LBA0926-h45). The 
first three objects all contain a DCO. These results con- 
trast strongly with the results in HOI and G09 in which 
the CHAf 036.3 absorption-line in the FUSE data is black 
or nearly black at line center in all cases. 

To understand our approach in using the CII line as 
a probe of the Lyman continuum, let us consider two 
limiting idealized cases. In the first, we assume a picket- 



fence model in which the far-UV source is surrounded by 
a population of clouds that are optically thick in both 
the absorption line and the Lyman continuum and that 
cover a fraction fc of the sky as seen by the source. In 
this case the residual relative intensity at the core of the 
line (/q) is just la/Icont — ^ ~ fc- Thus, a significant 
residual intensity in the line core implies a significant 
escape of Lyman continuum radiation. 

In the second case we assume that the far-UV source 
is surrounded by a uniform shell with unit covering 
factor. As discussed in HOI, in this case the optical 
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depth in the ISM just below the Lyman edge is related 
to the optical depth at the core of the absorption-line by: 

1) TLy/To = RZ-^\{agas/WO km s"!) 

where Zgas is the gas-phase carbon abundance scaled 
to the solar value of log{C/H) + 12 = 8.43 (Asplund ct 
al. 2009), we assume that Cll is the dominant ionization 
state in the neutral gas, agas is the velocity dispersion 
in the neutral gas, and R = 17 and 12 for the 1036.3 
and 1334.5 lines respectively (Morton 2003). Since the 
optical depth at the Lyman edge is over an order- of- 
magnitude larger than in the line core, in this idealized 
case a significant residual intensity in the core of the 
absorption lines is a necessary but not sufficient condition 
for the escape of ionizing radiation. 

To test these two models, we have used the set of Sill 
absorption-lines available in the COS data. Like Cll, 
Sill is the dominant ion of silicon in the neutral ISM. By 
using lines that all arise in the same ionic species, but 
which span a range in oscillator strengths, a comparison 
of the absorption profiles among these lines allows us to 
estimate the optical depth and covering factor of the gas 
(e.g. Heckman et al. 2000; Rupke et al. 2002; Quider 
et al. 2009; Erb et al. 2010). In the case of optically- 



thick clouds that do not fully cover the far-UV source, 
the line profiles of the different transitions should look 
the same. In the case of a uniform shell, the residual 
intensity will increase with decreasing oscillator strength 
/ (since tq oc A/). We will use the following transitions 
(with the respective values of logXf from Morton 2003): 
SiA111260.4 (-1.32), 1193.3 (-1.65), and 1190.4 (-1.95). ^ 
For each galaxy we simultaneously fit a single gaussian 
profile to the three different Sill lines, forcing the three 
lines to share a common centroid and width but allowing 
the normalization (depth) of each line to vary freely. In 
all four galaxies the results strongly favor the picket fence 
model with an implied covering factor of optically-thick 
clouds of about 40% for LBA0213+12 and LBA0926+45, 
50% for LBA0808+39, and 60% for LBA0921-h45. We 
show these fits for the case of LBA0921-I-45 in Figure 5.^ 
In principle, the sightlincs that avoid the clouds could 
be optically-thin to the metal lines but still be optically- 
thick to the Lyman continuum (equation 1). Thus, con- 

^ Unfortunately the Sill 1526 line (-2.18) is accessed in only 
one object and the 1304.4 line (-2.44) is badly blended with the 
SiIIIA1301.2,1303.3 stellar photosplieric features (Walborn et al. 
1995; Massa 1989) and the adjacent OIA1302.2 interstellar line). 

For the other four LBAs the data imply that the Sill lines are 
optically thick but fully cover the far-UV continuum source. 
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firmation that there is significant escaping ionizing radi- 
ation will require direct observations below the Lyman 
edge in these cases. 

In the mean time we can consider whether there is 
other evidence for the escape of ionizing radiation in 
these four LB As. In O09 we highlighted the relative 
weakness of the extinction-corrected Ha emission-line in 
most of the LBAs with DCOs, and speculated that a 
significant fraction of the ionizing radiation might be 
escaping. This information is listed in Table 2 where 
we compare the star-formation rates derived from the 
extinction-corrected Ha emission- line [SFRhci) to that 
obtained from the sum of the total-IR and far- UV lu- 
minosity {SFRiRuv)- It is noteworthy that three of the 
four objects above in which we have inferred an incom- 
plete covering of the far-UV source by neutral clouds have 
the highest ratios of SFRmuv / SFRhc- LBA0213+12 
(a ratio of 8.2), LBA0808+39 (3.0), and LBA0921+45 
(3.2). All three of these contain a DCO. In contrast, 
the mean ratio of SFRmuv / SFRhoc for the other five 
galaxies is 1.4. Thus, based on two independent lines of 
evidence, we suggest that there is a significant leakage of 
ionizing radiation from three of the eight LBAs in our 
sample. 

How do these results compare to what is seen at high 



redshift? There are four gravitationally-lensed high-z 
galaxies that have been observed with spectral resolution 
and signal-to-noisc similar to our sample. * Of these, the 
residual relative intensities of CIIA1334.5, 01 A1302.2, 
and Sill A 1260.4 (the strongest metal lines tracing the 
neutral ISM) are very small in MS 1512-cB58 (Pettini et 
al. 2002), The Cosmic Eye (Quider et al. 2010), and the 
8 O'clock Arc (Dessauges-Zavadsky et al. 2010), but 
large (^40%) in The Cosmic Horseshoe (Quider et al. 
2009). In fact, the last object appears to a good exam- 
ple of the type of picket fence situation seen in the DCOs 
in our sample. 

As discussed above, dust can also be a significant 
source of opacity to ionizing radiation in galaxies. Fol- 
lowing Shapley et al. (2006) and G09, we can distinguish 
between the absolute and relative escape fractions. The 
former is the actual fraction of ionizing photons that es- 

^ These strong interstellar metal absorption-lines from neutral 
gas exhibit significant residual intensity in the high signal-to-noise 
sta<;ked spectra of high-z star-forming galaxies (e.g. Shapley et 
al. 2003; Steidel et al. 2010). However, these data are obtained at 
relatively low spectral resolution (_R ~ lO''), so that the line profiles 
are not well-resolved in individual objects (and could miss gas with 
a high covering factor but low velocity dispersion). The stacking 
process itself could also smear out bla<;k, narrow components found 
in the spectra of individual galaxies. 
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Fig. 5.— The Sill 1190.4, 1193.3, and 1260.4 absorption-line 
profiles for LBA0921-I-45 all plotted on the same scale. In each 
panel the red line shows a fit of a single gaussian to the strongest 
component of the profile. This was a joint simultaneous fit to all 
three features, which were forced to have a common centroid and 
line width, but a free line depth. The similar fitted depth of the 
feature in all three cases shows that the data imply gas that is 
optically thick, but does not fully cover the continuum source (a 
picket-fence model). See the text for details. 

cape the galaxy, including the effect of dust. The latter 
neglects the effect of dust and is defined as the ratio 
of the fraction of escaping ionizing photons to escaping 
non-ionizing far-UV photons (e.g. it measures only the 
effect of the photoelectric opacity of the gas) . Our mea- 
surements discussed above provide information about the 
relative escape fraction. Given that the ratio of far-IR 
to far-UV fluxes for the galaxies with DCOs are of-order 
ten, the implied absolute escape fractions will be pro- 
portionately smaller. We list our estimates for both the 
relative and absolute escape fractions in Tables 1 (FUSE) 
and 2 (COS). 

4.3. Extreme Feedback 

As can be seen in Figure 4, the CIIA1334.5 absorption- 
line profiles in all the LBAs are blueshifted with respect 
to the galaxy systemic velocity, implying a large-scale 
outflow of gas. Further evidence for outflows is provided 
by the broad blue-asymmetric wings seen on the optical 
emission-line proflles (O09). 

The absorption-line proflles in Figure 4 show evidence 
for exceptionally high outflow speeds in the DCOs. To 
further investigate this, we turn to the Si III A1206.5 fea- 
ture. This is the strongest metal line that is accessed in 
all four DCO spectra, and arises in the ionized gas. As 
shown in Figure 6, outflowing gas is detected at extraor- 
dinarily high velocities in all four galaxies with DCOs. 



The flux-weighted line centroid is blueshifted by about 
700 km s~^ in these objects (Table 2). It is difficult 
to determine the maximum outflow velocity because of 
blending with the NIA1200 triplet, which lies ~1600 km 
s^^ blueward of the Silll line. Conservatively, the max- 
imum outflow speed seen in the DCOs reaches at least 
1500 km s"^ 

These velocities are much higher than in the local star- 
burst sample and in the LBAs without a DCO. This is 
shown in Figure 7 where we plot the outflow speeds in the 
ionized gas in local starburst and LEA samples vs. the 
star formation rate and galaxy mass (see Tables 1 and 
2). For galaxies with FUSE data we use the CIIIA977.0 
and/or NIIA1084.0 and/or CIIA1036.3 lines. 

The outflow velocities in the DCOs are also signifi- 
cantly larger than those typically seen in high-z galaxies. 
Steidel et al. (2010) find line centroids that are blueshifed 
on-average by 164 km and maximum outflow speeds 
that are typically ~ 800 km s""'^ for galaxies with star 
formation rates of ~ 10^ to IO^Mq year~^. 

Outflows from intensely star-forming galaxies are 
believed to be produced as gas clouds are accelerated by 
the ram pressure of a hot and fast wind driven by the 
collective thermal/kinetic energy supplied by supernovae 
and massive stellar winds (e.g. Heckman et al. 2000; 
Veilleux et al. 2005; Marcolini et al. 2005; Strickland & 
Heckman 2010), and/or by radiation pressure acting on 
dust (Murray et al. 2005; 2010). Can the unusually high 
velocities seen in the DCOs be explained in this way? 
We consider a simple idealized model of a cloud with a 
column density N accelerated outward from an initial 
radius tq by the ram pressure of a wind that carries 
momentum at a rate p into a solid angle fl. Then the 
terminal velocity of this cloud will be: 



2) Vte 



Here the momentum flux is in units of 10'^^ dynes, the 
initial radius is in units of 100 pc and the cloud column 
density is in units of 10^^ cm^^. The DCO star formation 
rates imply momentum fluxes of this order (Leitherer & 
Heckman 1995), while the HST images yield typical DCO 
radii of 100 pc (09). As discussed below in section 4.6, 
a column density of 10^^ cm~^ is a reasonable estimate 
for the ionized gas in these objects. 

The above neglects gravitational forces (cf. Wang 
1995; Martin 2005; Murray et al. 2010; Strickland & 
Heckman 2010). This is justifled here because the typ- 
ical sizes and radii of the DCOs imply virial velocities 
of only about 200 km s~^ (much smaller than the max- 
imal outflow speeds we observe). Note that Murray et 
al. (2010) conclude that outflows driven out of massive 
star clusters by radiation pressure will only have outflow 
speeds of-order the local escape velocity. This is not the 
case in the DCOs. 

These considerations imply that a galactic wind driven 
by an extreme starburst can in principle accelerate clouds 
to the high velocities that we see in the DCOs. The 
significantly higher velocities in these galaxies compared 
to typical starbursts can be naturally explained by the 
combination of a large number of massive stars confined 
within an unusually small radius. This leads to a very 
large wind ram pressure and/or radiation pressure at the 
launch point of the clouds. 
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Fig. 6. — Absorption-line profiles of the SiIIIA1206.5 feature. Contaminating foreground Milky Way and air-glow features are noted. 
Highly blucshifted absorbing material is detected up to at least -1500 km s^^ in the four cases with a Dominant Central Object (see text). 
Beyond this velocity there is possible contamination by the NIA1200.0 feature (as marked). 

we see in the DCOs (maximium outflow speeds of at least 
1500 km s^^) arc extremely rare in composite obscured- 
AGN/starburst systems and not seen at all in "pure" 
obscured AGN. They conclude that the AGN can drive 
very high-speed outflows, but these flows are confined to 
a small region very near the black hole (and are hence un- 
detectable in Type 2 AGN). They also conclude that the 
AGN plays little or no role in driving the large-scale out- 
fiows seen in the composite systems. It therefore seems 
unlikely that the high velocity outflows in the DCOs are 
driven by an obscured AGN. 

Tremonti et al. (2007) have detected outflows at veloc- 
ities ranging from 500 to 2000 km s~^ in massive post- 
starburst galaxies at z --^ 0.6. They speculate that such 
high velocities may require an AGN-driven outflow. It is 
also possible that the these objects are more massive ver- 
sions of the DCOs in this paper. This could be confirmed 
by HST imaging of their galaxies. 

4.5. Clues from, Lym.an Alpha 

The properties of the Lya profile are highly sensitive 
to the kinematics and distribution of the outfiowing gas 
and dust in star-forming galaxies (e.g. Stcidcl et al. 2010; 
Hansen & Oh 2006; Verhamme et al. 2006, 2008; Kornei 
et al. 2010). In Figure 8 we plot the Lya profiles for 



4.4. Feedback from an Active Galactic Nucleus? 

An alternative form of extreme feedback could be pro- 
vided by an AGN. As discussed in Overzier et al. (2009) 
the optical spectra of the LBAs with DCOs may imply a 
composite of both an intense starburst and an obscured 
(Type 2) AGN. The COS data show that the far-UV ra- 
diation in these objects is dominated by light produced 
by hot massive stars (Figures 2 and 3), which does not 
rule out the presence of a Type 2 AGN. 

Recent observations with XMM-Newton (Jia et al. 
2010) show that LBAs with such apparently composite 
optical spectra are also over-luminous in the 2-10 keV 
X-ray band by factors of roughly three to thirty com- 
pared to pure starbursts of the same far-IR luminosity. 
They conclude that it is likely that at least some of these 
objects do harbor an AGN, but that the overall bolo- 
metric luminosity of the galaxy is primarily due to an 
intense starburst. Could a Type 2 AGN be responsible 
for the unusually high outflow speeds and/or the appar- 
ently low covering factor of neutral gas due to a highly 
ionized state of the interstellar medium? 

Very high outflow speeds (many thousands of km s^^) 
are commonly observed in unobscured (Type 1) AGN 
(e.g. Rupke et al. 2005; Krug et al. 2010). However 
these same authors show that outflows similar to those 
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the eight LBAs with COS data, and highhght one very 
suggestive trend. 

In the five cases in which the inferred relative escape 
fraction is smah, there is very httle (if any) Lya emission 
blueward of the systemic velocity. Instead, the profiles 
show Lya emission redward of the systemic velocity and 
cither pure absorption or a mix of absorption and weak 
emission to the blue. Such profiles arc typical of high-z 
star forming galaxies (e.g. Shaplcy ct al. 2003; Stcidcl 
et al. 2010). In contrast, for the three objects in which 
we infer a high relative escape fraction (LBA0213+12, 
0808+39, and 0921+45) a significant fraction of the Lya 
emission is blucshifted with respect to the systemic ve- 
locity. To quantify this, we have measured the net equiv- 
alent widths (emission - absorption) blueward and red- 
ward of the systemic velocity in these eight LBAs. We 
define net emission (absorption) to have a positive (neg- 
ative) equivalent width. We then take the ratio of the 
red and blue equivalent widths (Reqw)- 

We list the results in Table 2. For the five galaxies 
with low inferred escape fractions Reqw is either nega- 
tive (redshifted emission and blueshifted absorption) or 
small and positive (strong redshifted emission and weak 
blucshifted net emission): LBA0055- 00 {Reqw = -0.21), 
0150+13 (-0.78), 0926+45 (0.14), 0938+54 (-0.17), and 
2103-07 (0.04). For the three galaxies with high inferred 
escape fractions Reqw is of-order one and positive (simi- 
lar amounts of red- and blue-shifted emission): 0213+12 
{Reqw = 0.36), 0808+39 (0.27), and 0921+45 (0.75). 

Profiles with significant blucshifted emission are un- 
common but not unknown in star- forming galaxies at 
2; ~ 2 - 3. Verhamme et al. (2008) discuss two such cases 
and interpret them as arising due to radiative transfer ef- 
fects in a static medium. This model can not be correct 
for the three DCOs in Figure 8 because the metal lines 



demonstrate that the gas is rapidly outflowing (Figures 

4 and 6). Erb et al. (2010) discuss the case of a high-z 
galaxy with both an outflow detected in the metal lines 
and a signiflcant amount of blueshifted Lya emission. 
They argue that in this case the optical depth to Lya 
photons is unusually low in the foreground outflowing 
material and that this is naturally related to the low 
covering factor of neutral gas they infer in this object 
(precisely the situation we see in our three galaxies). 

4.6. Where is the Dust? 

It is important to note that the observed UV spectral 
energy distributions (SEDs) of the DCOs are significantly 
redder than the intrinsic SED of a starburst (Ovcrzier ct 
al. 2011, and sec Table 2). For a normal starburst or 
SMC dust attenuation law (Calzetti 2001; Leitherer et 
al, 2002) these SEDs imply roughly two magnitudes of 
dust extinction along the line-of- sight to the DCO in 
the far-UV. This then raises a very interesting question. 
How can there be this much dust when the observations 
of interstellar absorption-lines from the neutral gas imply 
that such gas only covers a fraction of the far-UV source? 
The answer must be that this dust is associated with 
ionized gas that covers most or all of the DCO. 

This picture requires that the cores of the strongest 
absorption-lines due to the ionized gas should be nearly 
black (unlike the case of lines from the neutral phase). 
Examination of the profiles in Figure 6 shows that this 
is the case in LBA0808+39 (a relative residual intensity 
in the SiIIIA1206.5 line of 10%) , LBA0921+45 (7%), 
and LBA2103- 07 (<3%). In the case of LBA0213+12, 
the relative residual intensity in this line is high (~40%). 
However, the redshift of this galaxy is large enough to 
shift the CIIIA977.0 line into the COS G130L spectrum. 
This should be the strongest interstellar line due to ion- 
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Fig. 8. — Plots of the Lya feature. The location of the SiIIIA1206.5 feature is noted, as are the locations of foreground Milky Way 
features. In three cases (0213+12, 0808+39, and 0921+45) there is a significant amount of blueshifted Lya emission. These are the three 
objects that are candidates for a significant escape fraction of ionizing radiation. 

ized gas, and as such provides the most sensitive probe. 
Although the data arc noisy, the relative residual inten- 
sity in this line is small (10% or less). Thus, in this 
galaxy the Silll line is optically thin, while the CIII line 
is optically thick. 

We therefore conclude that the ionized gas has a near- 
unit covering factor in the DCOs. One plausible physi- 
cal model would be that the neutral gas represents the 
denser cores of clouds being accelerated by the wind, 
while the ionized gas represents a lower density halo of 
gas and dust surrounding the cloud core and which is be- 
ing ablated by the wind and photoionized by the intense 
radiation field (e.g. Marcolini et al. 2005). 

Can the DCO photoionize a column of gas with the 
required amount of dust? Elementary considerations 
of photoionization equilibrium imply that a source 
of radiation can photoionize a Stromgren slab with a 
column density given by: 

3) Nstrora = 1-2 X 10^3 U Cm'^ 

Where U is the dimensionless ionization parameter 
defined as the ratio of the density of ionizing photons to 
ek;c;trons in the slab. In order for a significant fraction 
of ionizing photons to escape, the actual gas column 



must be less than this value (to provide matter-bounded 

conditions) . The optical depth due to dust in the far-UV 
(1500 A) is related to the gas column density by: 

4) Ng^, = 3(5) X IQ-^^TFuvZ-^l cm'^ 

where the coefficients refer to an SMC (starburst) dust 
attenuation law (Leithcrer et al. 2002) and we assume 
that the dust-to-gas ratio is proportional to the metal- 
licity. The LBAs with DCOs all have gas-phase Oxygen 
abubdances close to solar (O09) and the line-of-sight red- 
dening implies Tpjjv ^ 2. Thus, to explain the linc-of- 
sight reddening towards the DCOs in a matter-bounded 
situation we need U > 7 x 10^'\ 

Alternatively, the optical depth at the Lyman edge 
due to neutral hydrogen in (mostly) ionized gas with a 
column density Nhii is given by : 

5) TLy = 8.6 X IQ-^'^NhiiU-^ 

For Nhii = 10^^ cm~^, we require U > 8 x 10"^ for 

TLy < 1- 

Is this plausible? Heckman, Armus, & Miley (1990) 
and Lehnert & Heckman (1996) showed that the ob- 
served radial gas density profiles seen in the optical 



12 



Heckman et al. 



emission-line nebulae in starburst winds are consistent 
with theoretical expectations for photoionized clouds 
exposed to the ram-pressure of the hot wind fluid 
(Chevalier & Clegg 1985). It is straightforward to show 
that in this case the ionization parameter in the clouds 
is given by: 

6) [/ = 9 X IO-3Q55/P35 

where P35 is the wind momentum flux in units of 10^^ 
dynes and Q55 is the starburst ionizing photon produc- 
tion rate (in units of lO'''"'' sec^^). 

The ionizing photons are produced by the most massive 
stars with lifetimes < 4 Myr, while the momentum in the 
starburst wind derives primarily from supernovae (whose 
progenitors have lifetimes of < 40 Myr). In equilibrium 
(e.g. a constant rate of star formation for longer than 
40 Myr), O55M35 - 0.5 and so [/ ~ 5 x lO-^. More 
generally, this ratio will be time dependent. For example, 
values oi U > 10~^ are associated with constant star 
formation for a duration less than 20 Myr and for an 
age < 4 Myr in an instantaneous burst. These ages are 
consistent with current constraints on the DCOs (O09 
and see Figure 3)^. 

We therefore conclude that it is feasible for the DCOs 
to photoionize a matter-bounded column of gas that is 
sufficient to contain enough dust to produce the observed 
reddening and the inferred attenuation of the far-UV con- 
tinuum and yet have a photoelectric opacity at the Ly- 
man edge that is less than one. 

5. SUMMARY AND IMPLICATIONS 

One of the major goals in cosmology is to understand 
the process by which the universe was reionized. The 
most plausible source for the ionizing photons is an early 
generation of hot massive stars in forming galaxies. The 
primary uncertainty at present is the unknown fraction 
of such ionizing photons that are able to escape from 
the galaxies and reach the neutral hydrogen in the inter- 
galactic medium. The column densities of gas associated 
with such intense star formation correspond to extremely 
high optical depths in the Lyman continuum. It seems 
nearly certain that some form of feedback from the mas- 
sive stars is required to produce the pathways through 
which the ionizing photons can escape. 

Our understanding of how this might occur would be 
greatly improved if we could find local examples of galax- 
ies from which significant ionizing radiation is escaping, 
and could then investigate in some detail the physical 
processes that enable this to occur. This will be most 
germane to the early universe if these local galaxies are 
good analogs to intensely star-forming galaxies at high- 
redshift. 

Accordingly, in this paper we used the Hubble Space 

Telescope Cosmic Origins Spectrograph (COS) and the 
Far-Ultraviolet Spectroscopic Explorer (FUSE) to inves- 
tigate the far-UV properties of a sample of eleven Lyman 
Break Analogs (LB As). These rare local objects are very 
similar to high-z Lyman Break galaxies in all their prop- 
erties studied to date (Heckman et al 2005; Hoopes et al 
2007; Basu-Zych et al.2007, 2009; Goncalves et al. 2010; 

^ We assume a standard Chabrier (2003) initial mass function 
and use the Starburst99 models (Leitherer et al. 1999) 



Overzier et al. 2009 - O09; Overzier et al. 2008, 2010, 
2011). We have compared these rare objects to a sample 
of fifteen typical UV-bright local starbursts. 

Following Heckman et al. 2001 (HOI) and Grimes et 
al. 2009 (G09) we have used the relative residual inten- 
sity in the cores of the strongest interstellar absorption- 
lines arising in the neutral interstellar medium (the 
C11A1036.5 and C11A1334.5 lines) as a probe of the ef- 
fects of the photoelectric opacity of the ISM on the Ly- 
man continuum. This approach is made possible by the 
high spectral resolution (R ~ 10**) and good signal-to- 
noise in our data. 

We find that there is a significant residual intensity (~ 
30 to 50%) in the Cll line core in four out of eleven of the 
LBAs but in none of the fifteen typical local starbursts. 
We show that in all four of these cases the residual inten- 
sity reflects an incomplete covering of the far-UV source 
by optically thick clouds (a picket fence model). Inde- 
pendent evidence for a signiflcant escape of ionizing radi- 
ation is provided by the anomalously weak Ha emission 
(relative to the sum of the far-UV and far-IR emission) in 
three of these four galaxies. It is therefore plausible that 
a significant fraction of the Lyman continuum is escaping 
in three of the LBAs in our sample. 

Assuming that the residual intensity in these lines does 
indicate the escape of ionizing radiation, then following 
Shapley et al. (2006) and G09, the relative escape frac- 
tion in these objects (defined as the ratio of fractions 
of escaping ionizing to non-ionizing far- UV photons) is 
~30 to 40%. After accounting for the effect of dust at- 
tenuation, the absolute escape fractions are ~ 3 to 10%. 

We have also examined the profiles of the Lya line in 
the COS data of the eight LBAs for additional clues. We 
find that there is strong emission both blueward and red- 
ward of the systemic velocity in the three objects in which 
a high escape fraction is likely. This strongly suggests an 
incomplete covering of the FUV source by neutral clouds 
(in agreement with the properties of the metal lines). 
The; other five cases show the profiles typically seen in 
high-z galaxies: Lya emission redward of the systemic 
velocity and either pure absorption or a mix of absorp- 
tion and weak emission blueward. 

In all three cases with high inferred relative escape frac- 
tions the HST images reveal that most of the far-UV 
emission originates in a single dominant central object 
(DCO). These are highly compact (radii ~ 10^ pc) and 
massive (one to a few xlO^ M0) young starbursts (O09). 
The interstellar absorption lines in all four of the LBAs 
with DCOs show extremely high outflow speeds: the line 
centroids are blueshifted by about 700 km s~^ and the 
gas reaches maximum outflow velocities of at least 1500 
km fi~^. These are signiflcantly higher than the veloc- 
ities seen in other low redshift starbursts or in typical 
high redshift galaxies. We have shown that these high 
speeds can be explained in principle by the acceleration 
of gas clouds in a starburst-driven wind. The exception- 
ally high velocities would be due to a large amount of 
momentum deposition in such an unusually small region, 
leading to an extremely high wind ram pressure and/or 
radiation pressure at the launch point for the clouds. 

The far-UV emission from the DCOs is significantly 
reddened and attenuated by dust. We have argued that 
this dust must reside in the ionized gas and have shown 
that the ionized gas fully covers the far-UV source (the 
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DCO). This situation is made possible by the intense 
ionizing radiation field in these galaxies (leading to sig- 
nificant column densities of dusty ionized gas that is op- 
tically thin in the Lyman continuum). 

Based on the above considerations, it is tempting to 
link the escape of ionizing radiation to the extreme feed- 
back associated with the DCOs. The combination of a 
high rate of energy deposition (kinetic energy and ioniz- 
ing radiation) within such a small volume may be the key 
to excavating channels through the interstellar medium 
and also intensely irradiating these channels so that they 
are fully ionized and translucent to ionizing radiation. 

It is interesting that only 19% (6/31) of the sample 
of LBAs imaged by O09 contain a DCO. Together with 
the present results this would suggest that large relative 
escape fractions are present in only 13% of LBAs. This 
is similar to the fractions seen in high- redshift galaxies 
(Shapley et al. 2006; Iwata et al. 2009). Could the 
presence or absence of DCOs in these galaxies explain 
these results as well? 

It is important to note that in the sample of LBAs 
studied by Overzier et al. (2009), the DCOs were prefer- 
entially found in the most massive galaxies (M* ^ 10"'^'^ 
to 10"'^^ Ms)- Galaxies this massive are exceedingly rare 
at 2: > 7, when the universe was reionized, and so the 
DCO host galaxies will not resemble the population of 
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low-mass systems believed to be responsible for produc- 
ing most of the ionizing radiation at these early epochs 
(e.g. Bouwens et al. 2009). 

We are instead suggesting that objects similar to the 
DCOs in this paper may exist in these early galaxies 
and facilitate the escape of ionizing radiation. A num- 
ber of theoretical papers over the past decade have in- 
vestigated the processes that regulate the escape of Ly- 
man continuum radiation from galaxies (e.g. Dove et al. 
2000; Clarke & Oey 2002; Fujita et al. 2003; Gnedin & 
Kravtsov 2008; Wise & Cen 2009; Razoumov fc Somer- 
Larson 2010; Yajima et al. 2010). These studies are at 
least qualitatively consistent with what we have observed 
in the DCOs. 

Finally, we note that O09 have shown that the masses 
and sizes of the DCOs are consistent with the regions of 
excess mass (cusps) seen in the centers of typical present- 
day elliptical galaxies (Hopkins et al. 2009). If 
so, then such galaxies may have all passed through an 
early stage during which they were sources of a signifi- 
cant amount of ionizing radiation that was able to reach 
the intergalactic medium. 
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